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I . INTRODUCTION 

This report presents the results of an analytical study 
of the MSFC Holographic Correlation Techniques Facility 
carried out under Contract No. NAS 8-323 62 by personnel of 
Montevallo Research Associates, Montevallo, Alabama during 
the period 12 January, 1977 to 12 June, 1977. Participating 
in the study were Dr. J. William Foreman, Jr. {Principal 
Investigator) and Mr. Joseph M. Cardone. 

The objects of the study were originally stated as 
follows : 

(1) To examine the existing MSFC experimental Holo- 
graphic Correlation Techniques Facility and to 
become familiar with existing theory pertaining 
to its operation. 

(21 To determine the theoretical effect, if any, on 
the correlation signal caused by insertion of a 
variable aperture into the object beam. 

C3) To examine theoretically the existence of a 
secondary correlation peak. 

(4) To determine theoretically the effect, if any, 
on the correlation signal caused by varying the 
surface roughness of the test object. 

(5) To explain theoretically why the object can be 
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translated in the object plane over distances large 
compared to a wavelength of light with negligible 
loss of correlation signal. 

It was anticipated from the beginning that it would 
probably not be possible to complete all five tasks during 
the present contract, and this turned out to be correct. 

A thorough theoretical understanding of the system, like 
a comprehensive experimental investigation of the system, 
will have to be developed as part of an ongoing study. 

Task (3) above was eliminated from the list when Dr. 

Robert B. Owen, our technical contact during the study, 
advised us that the secondary correlation peaks which had 
been observed earlier were apparently spurious. Tasks (1) 
and (5) have been completed, and a variation of Task (5) 
involving translation of the object along the optical axis 
of the system was begun. Time did not allow us to undertake 
Tasks (2) or (4) during the present study. 

It is suggested that the reader examine the Appendix 
before attempting to read the main body of this report, 
since many results which are used in developing the theory 
of the holographic correlation system are derived or summarized 
in the Appendix. 
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II. BASIC THEORY OF THE HOLOGRAPHIC 
CORRELATION SYSTEM 

The basic theory of the holographic correlation system 
is developed in a number of commonly available references . 
However, the usual treatment is somewhat abbreviated, and 
therefore difficult to follow in detail. We present here a 
treatment of the basic theory which is worked out in sufficient 
detail so that it can be readily followed. Our treatment 
follows the general approach used by Goodman in his dis- 
cussion of the Vander Lugt filter . 

Fig. 1 shows a schematic arrangement of a holographic 
correlation system. The operation of the system proceeds 
as follows : 

(1) An object mask with a certain desired transmission 
function is placed in the object plane, one focal 
length f in front of the transform lens. A 
photographic plate is placed in the hologram plane, 
which is one focal length f behind the transform 
lens. The photographic plate is exposed simul- 
taneously to the reference beam and the object 
beam. The exposure and development of the plate 
are controlled so that the transmission function 
of the developed plate is proportional to the 
total irradiance striking the plate during exposure. 
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Figure 1. Schematic diagram of a holographic correlation system. 
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(2) The developed plate, or hologram, is then replaced 
in the hologram plane , the reference beam is 
blocked off, and the hologram is illuminated with 
a new object mask placed in the object plane. 

The light throughput from the hologram passes 
through a correlation lens of focal length f 
' (identical to the focal length of the transform 
lens), and is examined in the correlation plane, 
which coincides with the back focal plane of the 
correlation lens . In general , four patches of 
light appear in the correlation plane, two of 
which overlap in the vicinity of the optical axis. 
One of the non-overlapping patches represents the 
cross-correlation of the original light amplitude 
distribution in the object plane (used to expose 
the plate) with the final light amplitude dis- 
tribution in the object plane (used to illumi- 
nate the hologram) . This is the term of interest 
here. Of course, if the hologram is illuminated 
with the same light amplitude distribution used 
to expose the photographic plate, the light patch 
represents the auto-correlation of this distri- 
bution with itself. 

Incidentally, the light distributions in the object 
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plane may not always be set up by passing a collimated 
light beam through a planar mask. They can also be set up 
by illuminating a planar reflector of some sort placed in 
the object plane in such a way that the reflected light is 
captured and processed by the transform lens. This is the 
way the MSFC holographic correlation setup operates at the 
present time. 

In order to set up the basic theory of the system, 
use will be made of the coordinate systems shown in 
Fig. 2. Coordinates in the object plane are denoted by 
1 *'»*>>' in the hologram plane by , and in the 

correlation plane by • Note that an inverted coor- 
dinate system has also been defined in the corre- 

lation plane. The reason for introducing this inverted coor- 
dinate system will be made clear later. 


Making the Hologram 

The first step in the theory is to describe the pro- 
perties of the hologram. Suppose the original distribution 
of light amplitude in the object plane is defined by the 
complex function 
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Figure 2. Coordinate systems for analysis of the holographic correlation system. 


8 


where and are the amplitude and phase, 

respectively, of the light vibration at a point in the object 
plane with coordinates Ike transform lens operates 

on the light leaving the object plane and, according to 
Eq. (A-17) , produces a light amplitude distribution 






in the hologram plane, where is the "lens operator" 
defined in Eq. (A-16) . 

A plane-^wave reference beam is also incident on the 
hologram plane. As indicated in Fig. 2, we will choose the 
direction of propagation of the reference beam to lie in the 

-plane and to make an angle 6 with the ■£- axis. A unit 

/v 

vector )k in the direction of propagation of the reference 
beam would have a negative component , as can be seen from 
Pig, 2, According to Eq„ (A-48) , the light amplitude dis- 
tribution in the hologram plane produced by the reference 


beam will be 




* * 


where t\ ft is the [constant) amplitude of the reference beam. 
Since we are dealing here with coherent laser light. 


the light amplitudes are additive, so that the resultant 
light amplitude distribution in the hologram plane will be 
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a, e * + ■]£ • 


The irradiance distribution in the hologram plane, I Ux, > 
is thus 


= 




ft. + 


(M 


x I i^l 




° \ 




_ ■ 

<7 — 


( 1 ) 


To make the hologram, a photographic plate with suitable 
sensitivity and spatial resolution is placed in the hologram 
plane, and its exposure and development are controlled 
in such a way that the transmission TT of the developed 
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plate is proportional to the irradiance X (*Zj during 

'exposure. We thus obtain a hologram with a transmission 
function 

jC (.OU, =• »C * (2) 

where £is a complex constant. 

i 

Illuminating the Hologram 

After the hologram has been made, it is replaced in 
the hologram plane, the reference beam is blocked off, and 
a new light amplitude distribution is established 

in the object plane. The transform lens operates on this 
light distribution, producing a light amplitude distribution 

in the hologram plane. When the light wave producing this 
distribution passes through the hologram, the resulting 
light amplitude distribution in the hologram plane will be 

£ ^ 5 (* 2 ., J (**, *§x) . 

Using Eqs. Cl) and (2) we then have 


it (* 2-5 = 


c 


U 


5 X 
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+ ^ e * 

-A. CTjU>^i 

+ ~ lz <$.0 J (*V e ^ ( 


Finally, the correlation lens operates on this light amplitude 
distribution and produces a new light amplitude distribution 
W ('Xj, *^. 3 ') i n t ^ ie correlation plane given by 

Wt-Kj,^ = -t- [&(Xt, <>»}]. (4 , 

Putting Eq. (3) into Eq. (4) gives 

w(x,,^ = X l3 [ 1 (**., 

(wr 

+ tzi [ Sl*.,4j0 

inr ) J 

• nr ^4^ $ 

+ s(^,>)e * ] 
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4“ 



5 tot* 



'X 



(5) 


For convenience, let us label the four terms in Eq. C5) as 
follows : 






( 6 ) 


^ (*», *&») * Tvfyt- t ~ t * t ’ ^ V* 1 ’ ^ v\| 


wj to*, = " ^L . ^ 3c^*-> ^ tot, *j-0 < 

(^4 ) 


(7) 

y ~r~ 


( 8 ) 


] 


nryyA^.6 

e * 

( 9 ) 


] 


Eq. C5 1 then becomes 


W^jAAj) s VJ_\ 4- £ 3 ) +- w s t^?^) ■*■ 

w f V (10) 
We will now investigate the individual terms in Eq.(10) in 

detail . 



First Term 
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From Eg. (6) , this term is 


But 









(6) 


Eg. (6) therefore becomes 


^(*3,^3) - - t z * ^ [& 

But from Eg. (A-37) we have 

= (wf 2 : (-**-&)■ 

Eg. (11) thus becomes 

* S ft . A (13) 


From here on, it will be convenient to work in terms of the 
inverted coordinate system in the correlation plane 

in order to avoid negative arguments. Since 
, we can write Eg. (13) as 




= £ 0. ■ a*) 

The light amplitude distribution in the correlation plane 
represented by this team is centered at the origin and is 
proportional to an inverted version of the light amplitude 
distribution £1X1,0 used to illuminate the hologram. 
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This term is of no interest in the present context. 

Second Term, Wz 

From Eg. (7) we have 

* ( 7 ) 

Using the convolution theorem stated in Eg. (A-32) we can write 


£;>*£ 




t 


But from Eg. (.A- 37) we have 




(15) 


^ ~ 'tt* \ °Cw. [ 


■= $>(**> *>*) ‘ 


( 16 ) 



15 


and from Eq. (A-43) we also have 


a ^l) 




*J-3 


(17) 


Putting the results of Eqs. (12), (16), and (17) into Eq. (15) 


yields 


t* \ 1 (**» ^ >) 4 (*^> v^l 

= (^V* ^s) * £*(-#*>-V») * #.(*»>**) • 

(18) 

We will now evaluate the convolution <j^(— #3 j..a^^4 ^ (X. 3 ^ 3 ) 
using the procedure outlined in Sec. B of the Appendix. 


^ / j \ 

Change variable names: t (-*>,-&) 

£ (*3> Vs) 

Fold £*(-f,-«l) i £(rS,-\) 

Translate <£&%) by amounts 
Kj and in the 5 and'YL 

directions , respectively ; x (sa) - 


^(-J.-'n) 
$. Ua) • 


£*(s-*»\vV s) . 
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Multiply the functions 
and integrate : 4. eo 

[[ £*(!-< vs-0 ** • 
— 00 

Our final result is thus 


+ *> 


J J (19) 

- *0 

Comparing Eg. (19) with Eq. (A-4a) we see that the integral 
in Eq. (19) is just the auto-correlation of with 

d*1iS6Xf * 

Butting Eq. (20) into Eq. (18) we then, have 


* ® £(x( *£3)] . (21) 

I 

Finally, putting Eq. (21) into Eq. (7) gives 


The light amplitude distribution in the correlation plane 
represented by this term is centered at the origin and is 
proportional to the convolution of i(<^) with the auto- 
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correlation function of * This term is not of direct 

interest in the present context. 

Third Term, 

From Eq. (8) we have 

. nrtf z / ov*. 8 

(M ) 3 

The terms within the square brackets in Eq. (8) represent 
the light amplitude distribution in the plane caused 

by the light wave leaving the t ) plane on its way to 

the correlation lens. It follows from the discussion in 
Sec. E of the Appendix that this light wave is travelling 
in the same general direction as a plane wave which would 
cause an amplitude distribution of the form g* \ 

in the plane. This direction, according to Sec. E 

of the Appendix, is along a unit vector which lies in 
the plane, makes an angle S with the 4- axis, and has 
a positive x^ component. This information is summarized 
in Fig* 3. Since this light wave is travelling at an angle 
g with the optical axis, we would expect the patch of light 
formed in the correlation plane by the correlation lens 
acting on this light wave to be centered off-axis in the 
correlation plane. This is the case, as the theory will show 
and as Fig, 3 indicates. 




Figure 3. Geometry of the light beams forming the various light 
patches in the correlation plane. 


i-* 

oo 
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Using the convolution theorem stated in Eq. (A-32) 
we have 

r : rrx a Av-> e 


(M\ 


{Vrt 


r "i r 2, b 

t j^e ^ 


4s^V&o]* *»»\ e *S& s ® 1 ] 


T~T ^ ^vi ^ v?. 1 1 

l/*v) 


* 


t*\ 


1 2.T& &va & "1 
e * 1 


= £(-Xs, -*&*)* ^ ~‘T^j • (23) 


Putting Eq, (A-46) into Eq. £23) gives 

* 0'% (*s> v) * n (*' , ,Q 4i J(j' + U^), (24) 
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Putting Eq. (24) into Eq. (8) gives 

- -=-£*- *<.**,■&)* 

vat) 

If we write the last equation in the form 


= . rr . & £ (^^4^6, ^.3) , 

• (H) J (25) 

then, in light of the discussion in Sec. C of the Appendix 
and Eq. (A- 12) in particular, we can interpret Eq. (25) as 
the convolution function 


£n. 

(U) 




centered at the location of the delta function £ (^3 ^ ® > ^ 3 ) 3 

which is at ^3 ® — 4 ^ S 0 . This term is sometimes of 

interest in optical filtering applications, but it is not 
of direct interest here. 


Fourth Term, 


Prom Eq. (9) we have 


ariw.6 




(W'l 

Proceeding as usual, we have 


( 9 ) 



p j. nru 
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&* £ 3f tt j\ ^ 


m * 6 


] 


* 4s [lWi,«j.ti] * <«[« ^ *r ] 

= ^**i*4-'lQ^ j ■* t~ 

* U [ 


Ai . 

s % (fajtys)* -^3j , 


where use has been made of Eg. CA-47) in the last step. 

Changing the arguments over to 4 ^ and we then have 

*»UXv> *(*», ^) C‘ * 

s ( ^ k 4 ) ^j^(— W 3 ■) “^ 3 ^ ^ £. (^3> ^ 3 ) ^ ^ 3 ) ' (26) 

Now from Eg. (20), we see that 

* i < 27 > 


Putting Eg. C27) into Eg. C26) gives 
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vb r i -) 

*3.3 jECXt^^ e * ^ J 

= (xtf [ #&»,<&) ® i i<$] * r (4-^4,^). 


(28) 


Putting Eg. (28) into Eg. (9) then gives 


* 7 ®£r ^3)^)1 4 S(o?3— 

(WO J (29) 
The light amplitude distribution in Eg. (29) may be inter- 
preted as the cross-correlation function 


~ y- 


centered at the location of the delta function 
Which is at ^"=4^$ j eft. This is, of course, 
of interest in holographic correlation systems . 


the term 


The term 



in the sguare brackets in 


Eg. (91 implies that the light wave leaving the ( ^ 2 ^ 2 .) 
plane on its way to the correlation lens, and which will 
eventually produce the light amplitude distribution in 
Eg. (291 in the correlation plane, is moving generally along 
a direction specified by a unit vector /X which lies in 
the -plane, makes an angle Q with the -&-axis, and has 
a negative x^-component . As shown in Fig. 3, this wave 
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(called the "correlation beam") leaves the hologram plane 
travelling generally parallel to the original reference beam. 
The correlation lens processes this wave to produce the 
"correlation spot" W 4 . centered at ^ « •£ 4 ^ & ? -s 0 . 

In order to separate the correlation spot 
from the two spots W\ Cfc*,., and W 3 . centered at 

the origin, it is necessary that the angle 0 at which the 
reference beam comes in with respect to the £ -axis be made 
sufficiently large. In the present MSFC holographic corre- 
lation system, the functions and are very 

small patches of light in the object plane, so that the 
spots W V > and 

are very small in size (less than one millimeter diameter) . 
Consequently, if cm, there would be plenty of separa- 

tion between the correlation spot and the two central spots 
which are of no interest. Since f is approximately 15 cm 
in the MSFC setup, the condition on 0 would be 


0? tw) ^ ^ ( \ cm) , 

which gives 

6 (jsf) -3,8 . 

Conversely, the angle 9 cannot be made too large, or 
else the correlation beam will miss the correlation lens 
and the system cannot operate properly. With reference to 
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Fig. 4 it will be seen that in order for the correlation 
beam to just strike the outer edge of the correlation lens, 
we must have 

^ - aw i } 

or 

6 $ W‘ ( • 

In the MSFC system, ^ LENS is about 10 cm, so we must require 
that 

0 ^ W ' ("||'') = s? -' 1 • 

In the present NSFC setup r the reference beam angle 
therefore should lie in the approximate range 

« 36 . 

Finally, it should be mentioned that the photodetector 
Which monitors the correlation spot in the correlation plane 
responds to the irradiance in the spot, not to the amplitude 
of the light vibration. The photodetector output current 
i will therefore be proportional to 

J [ I I Ax' 

A 



CeiWeLKTiftA 

flAMF 





Figure 4. Geometry to determine the maximum 
possible value for the reference 
beam angle 0 . 
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where A is the area of the photodetector aperture, under- 
stood to be centered at ^ Av*' 0 j = 0 . 
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III. TRANSLATION OF THE OBJECT 
IN THE OBJECT PLANE 


We now want to investigate the effect of translating 
the object distribution R in the object plane when the 

hologram is illuminated. The plate is exposed using the 
distribution , so that Eqs . (1) and (2) remain valid. 

Suppose now, however; that the hologram is illuminated with 
the function — W*) centered at ^ e. a b instead 

of the function JV centered at the origin. When 

the transform lens operates on the distribution 
the resulting light amplitude distribution in the hologram 
plane will be 

But using the shift theorem for the operator expressed 
in Eq. (A- 26 } r we have 

“i 




e(\% • 


The light amplitude distribution in the hologram plane thus 
becomes 


| ^ (Ma+b 










When the light wave producing this distribution passes 
through the hologram , the resulting light amplitude ‘ 
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distribution in the hologram plane will be given by Eq. 
(3) with replaced by 






The result will thus be 








5. ft. 


(H'i 


- VM*) I e 

fc 




2.TT 




0*)' 

The fourth term of Eq. (30) leads to the correlation spot, 
so the only term of interest to us in Eq. (30) is 


**(*vv) " v* 15 ^ iw*, v> 

(w) 


X t, 




. ( 31 ) 


(30) 
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When the correlation lens' operates on this term, the resulting 
light amplitude distribution in the correlation plane will be 


2 (*»,»*») 

J * (32) 


* e 


Use of the convolution theorem . (A-32)^j as usual then gives 



From Eq. (27) we have 

= $(***$») ® £rW>^)* (34) 

In order to handle the last factor in Eq. (33) , we can make 
use of an obvious generalization of Eq. (A-47) , 

&*[*> , (Wff (-**-♦=, (35) 

Comparing Eqs. (35) and C33) we see that 

C a Am\ & -V — 

4 


so -that ye have 



s £ (— ^3 & - & •) —^,3 — 



= ( ^ ■) £ (^3 — fat * 0*^ — k) 


(36) 
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Putting Eqs. (34) and (36) into Eq. (33) gives 

V) * ® J3r * S' W*-$**G - ^ 1 % -k) * 

CW (37) 

Eq. (37) states that the light amplitude distribution 

*$%) consists of the cross-correlation function 

Ag- 

centered at the coordinates a ~ ^ - Comparing 

this result with our previous result in Eq. (29) , we see 
that the only change has been a translation of the function 
tyjy by amounts a and b in the Xg and y^ directions, respectively. 
Thus, translation of the object distribution in the object 
plane merely translates the cross-correlation function in 
the correlation plane. Movement of the photodetector to 
the new location of the cross-correlation function should 
then produce the same output current. 

The effect of translating the object distribution in 
the object plane is summarized in Fig. 5 for convenient 
reference . 



o&oect 



illuminated with the object 
distribution 



(b) Situation when the hologram is 
illuminated with the translated 
object distribution 


Figure 5. Effect of translating the object distribution 
ft in the object plane. 
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IV. TRANSLATION OF THE OBJECT 
ALONG THE OPTICAL AXIS 

As the final task in the present study, we started 
to look at the effect of translating the object distribution 
along the optical axis of the system (i.e., along the 
z-axis in Fig. 2) . We were not able to complete this task, 
but we report here the partial results we have obtained. 

Here again we suppose that the photographic plate is 
exposed using the distribution , so that Eqs . (1) 

and (2) are valid. Suppose, however, that when the hologram 
is illuminated with the distribution ' t * ie object 

plane is at a distance f+A£from the hologram plane, instead 
of the usual distance f. Then, according to Eqs. (A-13) and 
(A-16 ) t the light amplitude distribution in the hologram 
plane will be 



in Eq. (A-13) as usual. When the hologram is illuminated 
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with this light amplitude distribution, the resulting light 
amplitude distribution in the hologram plane will be given 
by Eq. {3) with replaced by 


-* 1 - 1 **} i(M0 . 


The result is 


Cn 


s(* 4)1 ^e '*+ 




WP 



so let us concentrate on this term; 


£t»,v>2»s „i 


(38) 
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When the correlation lens operates on this term, the resulting 
light amplitude distribution in the correlation plane will be 


*%(*>. <»[&<*•. *) S.Kfr) 

r 


x e 


_ i ?. few. 6 

9 K 


] 


Application of Eqs. (A--32) , (A-37) , (A-43), and (A-47) then gives 


= £ h * * 


(H) 


3 *(*.*»> 


(39) 


We have to pause now and calculate the third factor in Eq. (39), 
By definition of the operator ££ given in Eq. (A-16) , 


we have 


^ J J • 

(40) 


- 00 


Now since ^ a 32£ , we have 




* e 
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= e 


-qg- «^\) 


4 


Thus Eq. (40) becomes 








“ *6 


4“06 

B \\ e%F^'- + V + <M* ■•• vfrO] 

— 06 


. (41) 


Completing the square in the argument of the exponential 
term in Eq, (41) gives 


TTA^ 

e ' t W r 



+- 

£& 



* e 


^ ^ C V -v *3) + (*^- 


K*- 


x e 


T^“ 



Eq, (41) can thus be written as 
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* * Ax, A 


« e 


— oo 


(42) 


The integral in Eq. (42) is just a constant, the value of 
which is obtained in Sec. F of the Appendix, Eq. (A-65) . 
We thus have 

a .it #-*%.>■ 






&rb 




(43) 


Putting Eq. (43) into Eq. (39) gives 




^ 4 As 






= £ ( ~*z *-«£) * *i) * eis ( * 3 ' + ^ ) 

X-t Ai 


* JC*3— ■j.,') . (44) 



37 


Use of Eq. (27) then gives 




£ ■>» 

42r 





* $(*,'- •f fi*YS ^ ^ »'o. (45) 


The light amplitude distribution in Eq. (45) may be interpreted as the 
function 


5 ^ 
to A* 





centered at = f sin 0 , y^ ~ 0. This function represents the 

usual cross-correlation function, 
convolved with the spherical phase factor 


e Vi 


Unfortunately, this is as far as we have been able to develop the 
theory within the time frame of the present contract. 
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V. CONCLUSIONS AND RECOMMENDATIONS 

We regret that the study of the effect of translating 
the object along the optical axis reported in Sec. IV could 
not be completed during the present contract. In practical 
holographic correlation systems used to examine three-dimen- 
sional objects, it is certainly possible for the objects to 
change size, or to undergo small translations, along the 
optical axis. We would recommend that this study be brought 
to completion eventually. 

An observation which strikes us as being very impor- 
tant is the fact that while the theory of the holographic 
correlation system is always developed on the assumption that 
.the object plane, the hologram plane, and the correlation 
plane are all parallel to each other, in practical setups 
this is usually not the case. The effect of having these 
three planes slanted with respect to each other has not been 
analyzed in the literature, as far as we have been able to 
determine. We feel that this is an important study which 
ought to be given highest priority. 
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APPENDIX 

The purpose of this appendix is to summarize some impor- 
tant results which will be needed in developing the theory 
of the holographic correlation system. 

A. The Fourier Transform 


Let be a complex-valued function of two indepen- 
dent real variables, x and y. (NOTE: In this report, complex- 

valued functions will be indicated by an underline.) The 
two-dimensional Fourier transform of is defined as 

where and £1^ are spatial frequencies with dimensions of 

and •— , respectively, and The inverse Fourier 

transform of T is defined as 


— CO 


B. The Convolution and Correlation Integrals 

The convolution integral of two functions and 

is denoted by ' and is defined as 


follows : 
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-Y-oo 

= U k\- (A-3) 

- oO 


It will be observed that the result of the convolution operation 
is a function of the original variables x and y. This function 
is called the convolution function of ^ and 
5 an< ^ in Eq. (A-3) are just dummy variables of integra- 
tion. 

The convolution operation is' both commutative and associative. 
That is. 


<L. . The variables 




and 


Since the convolution operation is both commutative and 
associative, it follows that the individual operations in 
a string of convolutions can be performed in any order whatever. 

In order to set up the convolution integral of two 
functions and (.*>$•)' the following steps are taken: 

(1) Change the names of the variables from (x,y) to 

<5-n> • The functions then become and 3 ^,%) . 

(2) Fold about the origin to get . 

(3) Translate by amounts x and y in the 

and ^ directions, respectively, to get ? , 
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(NOTE: Since convolution is commutative, steps 

(2) and (3) could equally well have been applied 
to the function • The final results would 

be identical either way.) 

(4) Multiply the folded and translated function 

*>v0 by the other function t and. 

integrate over all space to get the final result. 


\ — 

= [ \ •yx) *5 ^ • 


— o o 


The cross^correlation integral of two functions 
and fjg&ytfri is denoted by ® , and is defined as 


follows : 


oo 


= U is i\ • (a-4) 

— od 

Notice that the resulting function, called the cross-corre- 
lation function of^ and , is a function of the original 
variables x and y. 

The cross-correlation operation is not commutative. 
That is, 


However, the translation operation which has been applied to 
& 

the function in Eq. Ch-4) can be applied instead to the 
function ^in tfle opposite sense without changing the value 
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of the integral . That is , 

+ o© 

® %CVj) = \\ ^.*(1,^ ^ : (3+*,n+'3-) (a-5) 

— oO 

(NOTE : If the functions ^CXjA/p and 4 l M) are identical, 

then the operation 


4-o* 


(A-4a) 


— <3* 


is called the auto-correlation function of ^ with itself.) 

The cross-correlation integral of two functions 
and can he set up by use of the following steps : 

(1) Change the names of the variables from (x,y) to 
obtaining the functions j^.($', /T 0 and U) • 

(2) Take the complex conjugate of the first function 

jjf. 

named in the operation, obtaining . 

(3) Translate the function ^ \5;‘ T l) by amounts x and y 
in the ^ and ^ directions, respectively, to get 

£*(?-*, vVb 

(4) Multiply the conjugated and translated function 

by the second function ’ an< ^ 

integrate over all space to obtain the final result. 




^.oO 


— o© 
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Notice that convolution and cross-correlation differ 
in two respects, folding (done in convolution but not in 
cross-correlation) and complex conjugation (done in cross- 
correlation but not in convolution) . 

C. The Dirac Delta Function 


The well known Dirac delta function has several impor- 
tant properties which we wish to summarize here for convenient 
reference. For purposes of this report we are interested in 
the two-dimensional version of the delta function. In the 
equations which follow, a and b represent arbitrary real 
constants . 


Basic Definition 




00 if x=y=0 
q otherwise 




4- oa 


J $ (ft-OCjb-Aj.') 


~ 1 , (A-6) 


— 00 

Sifting Property 

4 - ot> 




(A- 7) 


w po 


where j|Mx r y) is an arbitrary function of x and y. 
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Scaling Property 

£ (a* , 1> = — ■ . (A-8) 

|«.W\ 

NOTE: It follows from the scaling property expressed 

in Eq. (A- 8) that 

S Ox-*, - i 

* 1 _ 



so that the delta function is a symmetrical 
function of its arguments, and Eqs. (A-6) and 
(A-7) can equally well be written 

+ o* 




1 (A- 9) 




— cO 

Integral Representation 


v C&, b) . < a - 10 > 


+«■«> 

[[ 




^ (ft 4 


+eO 


i i ff -4 2nr (*** + -‘y^ 

= J \ e AxA 




3- 


= S • <*-n> 
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Convolution of a Given ^Function with a- Delta Function 

Suppose a function is convolved with a delta 

function located at . Following the four 

steps listed earlier for setting up a convolution integral, 
we have: 

(1) Functions become and 

(2) Either function can be folded and translated; 

we choose to fold' and translate • Folding 

produces 

(3) Translation produces ^ (0C- f , u ^~ . 

(4) Multiplication of the functions and integration 
over all space gives 

= J | £(*-?, <j.-V) i (S-a, V*>) H- 

- oO 

Using the sifting property of the delta function expressed 
in Eq. (A-lO) ,"-we now have 


5 « *» 

a g (0C-<k 7 aj.-W') . (A-12) 



46 


The result, according to Eq. (A-12) , is a translation of the 
function by an amount a in the x-direction and an 

amount b in the y-direction. Stated another way, we can say 
that the function has been translated to the location of the 
delta function, x=a and y=b. 

It is difficult to indicate graphically the effect of 
this operation in two dimensions, but a one-dimensional picture 
is feasible. The diagram shown below indicates graphically 
the result of convolving a real function with the delta 

function £(%— a} located at x=a. 



D. The Transformation Properties of a Lens 

Suppose a light wave advancing along the -^-.axis gives 
rise to a light amplitude distribution in an object 

plane at a distance in front of a lens of focal length f, 
as shown below. When the light wave is acted upon by the lens. 
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In the special case where the object plane coincides 
with the front focal plane of the lens, so that&=£, the 
factor in the exponential in front of the integral 

in Eq. (A-13) vanishes, so that the exponential becomes equal 
to unity, and Eq. (A-13) becomes 


4* od 


= -p- 1 e 




‘ (A-14) 


— cxO 
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Now the factor — s-^ = ^ in front of the integral in 


Eq. (A-14) represents a constant phase shift of-3£ , which 
is of no practical importance. We can thus drop this factor 
and write 


Vr - j- *+ 

e Ax.A-w, . (A-i5) 

The integral in Eq. (A-15) can, of course, be written 
in the form of a Fourier transform by defining spatial 
frequencies -P* = V« and % = . However, it is our 

opinion that when this is done and the theory of the holo- 
graphic correlation system is developed using Fourier trans- 
form notation, considerable confusion can result about the 
arguments of the functions involved, and constant factors of 
M tend to be dropped indiscriminantly . 
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The theory of the holographic correlation system can also 
be carried through using the integral operation in Eq. (A-15) 
directly. This is the approach used by Smith 3 . However, this 
procedure becomes very cumbersome, and in the process of using 
this approach one essentially derives anew many of the commonly 


known properties of the Fourier transform. 

Because of these problems, we will use a somewhat 

t 

erent approach here. Let us define an integral operator 
f, based on Eq. (A-151, as follows; 

v-tf - -T[ *«■»> 

(A-I?) 


The symbol " ££ " stands for "lens operator," and the subscript 
"■12" indicates that this operator transforms a light amplitude 
distribution^*,,^,) in a plane at a distance f in 

front of the lens into a light amplitude distribution 
in a plane at a distance f behind the lens. From 

Eqs, (A- 1 5 ) and CA-16). , we can write the distribution 
in the back focal plane of the lens as 


(A- 17) 

We will now proceed to establish those properties of the 
operator which will be needed in developing the theory 
of the optical correlation system. These properties, of 
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course, are analogous to the corresponding properties of the 

standard Fourier transform. 

-I 

The inverse Operator, ^Y3» 

-l 

To establish the form of the inverse operator , 

we will first convert Eq. (A-16) over into a standard Fourier 
transform. If, as indicated earlier, we define spatial 
frequencies 



* 


* 


0^ (A— 18) 


, (A- 19 ) 


then Eq. (A-16) can be written in the form of a standard 
Fourier transform. 


r T r + rt -j **(*&< 

*n.[***\ -- ***>*>- % KWj) = ] ^ 


— et 


(A-20) 


From Eq. (A-2) , the inverse operator would be 

-t, 


4 [V*-*9 - - 3- [l = S Jit . 


—ot 


(A— 21) 
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Now from Eqs. (A-18) and (A-19) we have 


= k { *4) 

= k ( ’ 


and we also observe from Eq. (A-20) that an< ^ X 
are just different names for the same function. It then 
follows that Eq. (A-21) can be rewritten in the form 






-t£> 


MfrW 


M- ) 


4 00 


on 


Jj#*^ 


i-fr (v^+>-p 


•z. 

— po 


A'Xa.A^ . (A-22) 


-t 

Eq. (A-22). establishes the form of the inverse operator t ix . 
Shift Theorem 


Suppose that 


r -1 + ft - j- 

t, x l e 


— «o 


(A-23) 
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Then if cfu. operates on the translated function 


we have 




— oO 

If we introduce new variables 

C<= 

f = ty-b, 




A^(| * 

(A-24) 


then we have 


K v = a- 

*&*■= f + *> 

A*, = M 

H 1 * 7 


and Eq. (A-24) becomes 




— . pd 


Act A|? 


= e 










— ** 


Act. Ap 

(A-25) 


Comparing Eqs. (A-25) and (A-23) , we see that 
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— {?© 


7 


since the name of the dummy variable of integration is 
irrelevant. Eg. (A-25) therefore becomes 




~H4 




. (A-26 ) 


Convolution Theorem 


Suppose that ^12- [ and 

= 5 [fa) * Then consider what happens when operates 
on the product function £tM>) • By definition 

of the &_ operator in Eq. (A-16) , we have 

+ r °f ~i 




But from Eq. (A-22) we can write 


(A-27 ) 


-i + ” s 




-«*a 


^ei ^6 

(A-28) ' 


where we have written the dummy variables of integration 
as oOj ^ instead of to avoid confusion with the variables 

Xz., 1 ^. 2 . in Eq. (A-27) • Putting Eq. (A-28) into Eq. (A-27) 
we have 



+ C*4 
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~cA 


(M) 1 


-SKfie A=a« 


- > ff- ( 1t,1(i ' , ») 

^ ^ A.*! A^., . 


Interchanging the order of integration gives 

+ et> 

Js£^)Ac<Ap 


-c* 






~ 'hf W [''‘ m (^2-°0 + *^l 




-(A 


(A-29) 


Comparing the final integration in Eq. (A-29) with Eq. (A-16) 
we see that 


+ <* 




- + D' 


Ax 


-</> 




Putting this result into Eq. (A-29) then yields 

+ 06 

j ^-f) U^. 


-A 


(A— 30) 
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Comparing Eq. (A-30) with Eq. (A-3) we observe that 


= (a ' 31) 

V«t) 

where the interchange in order is legal because the con- 
volution operation is commutative. But, by our original 


assumption. 


^ ° S \-2 1 . )\ 




Eq. (A-31) can therefore be written in the form 


^\2. \ 4> ^‘1 fy ) J2r L ^<2. 1 £ ^ h ^ S \ * ' 


(A-32) 


Double- Trans form Theorem 


What is 


Suppose again that ^ = ^*1,^0 • 

[ 4. I • 
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From Eq. (A- 16) we have 


-too 




(kX.A 


-33 ) 


■*iU^«- i ^‘ w ' , "‘ , W. 

“• «o (A 34 ) 


Putting Eq. (A-34) into Eq. (A-33) then gives 


+00 +oC 


— 06 




X e 


Interchanging the order of integration gives 


> 1 


4 [u*<»wa] = 5 \ *»■>*> I n e ’^ 


— mC> 


X 



Ax j Aj^ j * ( A 3 5 ) 


Now upon comparing the inner integral in Eq. (A-35) with 
Eq. (A-ll) we see that 



-oo 


fi 


eo 


f ^ 3 ) + ^-2 (ty + 





51 (^0 S (*l+*3 7 $i+$. 3 ) , (A-36) 


where the second step follows from the scaling property 
of the delta function expressed in Eq. (A-8) . Putting Eq. 
(A-36) into Eq. (A-35) we obtain 

-tea 

[^2. [ \\ $&***b') £ (*^3 > ^ 3 ) M H' * 

— oft 

Using the sifting property of the delta function in Eq. 
(A-10) , we finally obtain 

\t\v. [^* 1 ,^]] « (A-37) 

We will also have need for the double transform in 
a somewhat different case. In particular, if 


» oC\2 > 
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we need to 


evaluate 4 [ 4 ^ 0 ] . To begin with, from 


Eq. (A-16) we have 


V***' = [I £ ( M') « ” A*,Aj, ■ (A-38, 


It follows from Eq. (A-38) that 


# + n « 

. (A-39) 


" ©O 

From Eq. (A-16) we then have 

■f 


^ r * 1 ff .*,* \ ~fr 

* IU l ^ e 

Putting Eq. (.A-39) into Eq. (A-40) gives 


(A-40) 


U *‘‘M • H • * ” \\ *x,o 




x e 




Interchanging the order of integration gives 


-\\ *v,,> n 
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X 



(A-41) 


By use of Eqs. (A-ll) and (A-8) we then have 

•h c# 


!i 


— 


^ c$ 3 -9.0] 




S T-^r’ ^1 

(jX£) S ^ S -^0 * 


Putting Eq T (A-42) into Eg. (A-41) we obtain 

06 


— •a 


Finally, use of Eg. (A-7) gives 


Going back to the definition of , we 

write this result in the form 


(A-42) 




can then 


(A-43) 
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Effect on Complex Exponentials 


2-TTC.Xz 


Consider the operation 


where 


c is a real constant. From Eg. (A-16) we have 


'in ^ 


• arcXx. 

e*~ 


+•* 


\-\w 


ac *' -i m <**+*»«) 




~#6 


\ \ ^ ^ c ) 4- 

J i (Kx^m , • 


'C* 


Comparing Eg. (A-44) with Eg. (A-ll) we see that 


V 

(At- 44) 


II 


'cO 






«• 

- (^n £ C •> ■“^'3^ » (A-45) 


where the last step follows from Eg. (A-8) . Putting Eg. 
(A-45) into Eg. (A-44) gives 


^ = Cu)Vv^=, 


(A— 46 ) 
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; z-n-c^x 

^ J* * 

To obtain the result for cKzj operating on C 
it is only necessary to change the sign of c on both sides 
of Eg. (A-46) . We then have 


i z-rrcXx -u 

^ 1 ™ S (-OC3 - , -^0 * (A-47) 


E. Some Properties of Plane Waves 


Consider a plane wave propagating in a direction 

A, A 

defxned by a unit vector a. - We will take ik to lie in the 

A 

JC&- plane, as shown below, and we will suppose that AX makes 
an angle 0 with the£-axis and that the x-coxnponent of IX 



the light amplitude distribution in the xy-pl 


ane 


due to the incident plane wave moving in the direction JtX. . 
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The equation of a general plane wave is 

rue *-*" 1 

where is a constant amplitude and 



s 

Wk. 


where A is the wavelength of the light (assumed mono- 
chromatic) . From the figure above it will be seen that 


IK ~ £vrv 6,0, C<2-0) 

"ft = 0*? 0 , 

so that for any point (x,y) in the xy-plane we have 

IK-lC * 

A ^ ^ 

The light amplitude distribution in the xy-plane due to 
the plane wave in question is thus 

- £ 

^ . (A-48) 


Conversely, if a plane wave sets up a light amplitude 

distribution in the xy-plane given by Eq. (A-48) , 

, . A 

xt follows that the unit vector aK which describes the 
direction of propagation of the plane wave makes an angle 0 
with the -2-axis, lies in the plane, and has a negative 
x- component. 
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Following through the same reasoning for the case 
where lies in the plane, makes an angle 0 with the 
■2— axis, but has a positive x-component, one finds the light 
amplitude distribution set up in the xy-plane to be 

an rta#* e 

= n 0 e * # (a— 49) 


Conversely, if a plane wave sets up a light amplitude 

di s t r ibut i on in the xy-plane given by Eq. (A-49), 

~ A 

it follows that the unit vector IK describing the direction 
of propagation of the wave makes an angle 6 with the 
• 2 >— axis, lies in the plane, and has a positive x-compo- 
nent. 

By way of generalization of these ideas, Smith ^ 
shows that if a light wave sets up a light amplitude 
distribution 




* * 


(A-50 ) 


in the xy-plane, where is a general complex function, 

then the light wave causing this amplitude distribution 
is travelling generally in the same direction as the plane 
wave which caused the distribution in Eq. (A-48) , except 
that the wavefronts in the case of Eq. (A-50) are not planes, 
but rather curved surfaces of some type. Similarly, if a 
light wave sets up a light amplitude distribution 
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* 5 ! 


in the xy-plane, then this wave is travelling generally 
in the same direction as the plane wave which caused the 
distribution in Eq. (A-49) , but the wavefronts are curved 
rather than planar. 


F. Evaluation of the Integral in Eq. (42) 


Let us denote the integral in Eq. (42) by $ ’ 


+©£ 


il 




* 3 ) + + "] 


n 


** ** CAr-51) 


Let us also introduce a new constant 


a = 

^ n v 

and new variables of integration 


(A-52) 


e( = fc. 




(A-53) 


f * *&*+ • <a - 54) 

Since and $3 are constants as far as the integral in Eq. 
(A-51) is concerned, we have from Eqs. (A-53) and (A-54) that 
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Aot ss. 

= H* 


Furthermore, the limits of integration in terms of and 
are the same as the limits of integration in terms of 
and . Eg. (A-51) therefore becomes 


i . *fj uif 


4-o0 




-06 
•V od 


r + ** 

) J Cn.$,(A-f ) iMif - j. £avn 




+ <* 


* ii ***" ^*"1 


-“p& 


+ <* 


J[l** +- C#-<3,«< ^ (A— 55) 
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Now let us write, for convenience. 




(A-56) 


where 


+ «* 


^ — /W> a^U ^ 57 ^ 


“p* 


° fL ‘ ^ 


We then have 


tO» -V o<5 


Now consider the integrals 

•V a 6 


: v = ^ 


too 

^ CguQlcL*' 


+ «* 


> JUL 


s 2. 
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Z JL 

Introduce a new variable <*« • 
and we have 


Then 

4 <r 




r +eo 

I, = 


1 ^ 

X* = 


f +p6 

we have^ 



+ o* 
rs 


+ oft 

1 C<52.^ - 

A 


Thus 

+ ** + 04 

s j C ^.et = — |* Avw <3^o( 

— ©o 


= 


-oft 


TT 


Putting Eq. (A-60) into Eq. (A-59) we have 

+ 00 

,r . ..*■ nr 

Gif 


» W 

■» v = 


IT 

zq 


e 1- . pet 

r <i a* j 


— oft 


^ + *o + oo 

| i J ^ j * 


(A-60) 
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Eq. (A-60) then gives 

«i = 


IT 


l 




N 


— \ = 
2I * j 


O. 


(A-61) 


Going back to Eq. (A-58) , we have in a similar manner 

+ 4* ed +• od 

V*. - | ^ c ^- ^ 4- Avrv C<ft. <*oF * 

—p 6 


-pO ~ p & 

Invoking Eq. (A-60) gives 


4* o* 




•S 




TT, 

2Q 



( F 


T 


1 • ?• 


(A- 62) 


Putting Eqs . (A-61) and (A-62) into Eq. (A-56) gives 


i 



ir 
<* * 


(A- 63) 
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The result in Eq. (A-63) assumes that$>o since otherwise 

V 2- 

the change of variable would not be valid. From 

Eq. (A-52) , in turn, this requires that ^>0 . 

If , on the other hand , 6 , then <*<0 . Let us then 


define a positive real number such that 




(A- 64) 


We would then have 


ii 




-<A 


which following the same steps as before, would lead to 
the result 


% * * 

But in view of Eq. (A-64) , this becomes 


% = -*X. 

— ( n 

which is our previous result, 
is valid whether <^>5 or ^<0 
result 


In other words, Eq. (A-63) 

We thus have as our general 


~ 0 $ 


IT 


-if 


^ firdr£\ 


=■ — a, ^ 


Now — •£ is a constant phase factor of no interest, 

so we will drop it and write 


* = M 


(A-65) 



70 


REFERENCES 


1. J. W. Goodman, Introduction to Fourier Optics , McGraw- 
Hill Book Co,, New York (1968), pp. 171-177. 

2. A. Vander Lugt, "Signal Detection by Complex Spatial 
Filtering," IEEE Transactions on Information Theory, 
Vol. 10 (April, 1964), pp. 139-145. 

3. H. M. Smith, Principles of Holography , Second Edition, 
John Wiley and Sons, New York (1969), pp. 242-246. 

4. Reference 3, pp. 19-20. See especially Fig. 2.8, 
p. 20. 

5. See, for example. Standard Mathematical Tables , 19th 
Edition, Chemical Rubber Co., Cleveland (1971), 

p. 447, #636. 



